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Abstract: In this paper, we extend our previous analysis concerning the formation of 
inhomogeneous condensates in strongly-coupled fermion effective field theories on curved 
spaces and include the case of conifold geometries that represent the simplest tractable case 
of manifolds with curvature singularities. In the set-up considered here, by keeping the 
genuine thermodynamical temperature constant, we may single out the role that curvature 
effects play on the breaking/restoration of chiral symmetry and on the appearance of 
inhomogeneous phases. The first goal of this paper is to construct a general expression of 
the finite temperature effective action for inhomogeneous condensates in the case of four- 
fermion effective field theories on conifold geometries with generic Riemannian smooth base 
(generalised cones). The other goal is to implement numerically the above formal results 
and construct self-consistent solutions for the condensate. We explicitly show that the 
condensate assumes a kink-like profile, vanishing at the singularity that is surrounded by 
a bubble of restored chiral symmetry phase. 
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1. Introduction 

Current understanding of the phase structure of strongly interacting field theories is still 
far from being complete, despite the fact that this problem has been at the center of 
active interest for many years. The case of QCD is of particular physical relevance and 
many ongoing efforts are directed at mapping the geography of the phase diagram of 
the theory at finite temperature and density, understand the nature of the transitions 
between the different phases, locate the critical points, describe the morphology of the 
ground state. Impediments to arrive at the desired complete characterization of the phase 
structure of QCD are mainly due to the notorious difficulties in performing ab initio lattice 
computations under the influence of generic external conditions and, amongst the various 
approaches, the use of effective models is one of the most common playgrounds where the 
phase structure can be discussed (see Refs. [jl], ^, ||, [|, || for review). Within these effective 
models, the Nambu-Jona Lasinio (with its variants) is, possibly, the field theoretical set-up 
that received greatest attention. 

One particular aspect that has been at the center of recent discussions is the identifica- 
tion of inhomogeneous phases. In fact, the appearance of modulations in strongly-coupled 
fermionic systems has led to several new insights concerning inhomogeneous phases in QCD. 
More precisely, several works indicated that, in the chiral limit, inhomogeneous phases are 
energetically preferred for relatively large values of the chemical potential, and the preferred 
inhomogeneous ground state in the vicinity of the chiral critical point assumes a crystalline 
structure similarly to superconductors. The problem has been analyzed in different models 
and seems to be quite generic (see, for example, ||]). The inclusion of gauge degrees of 
freedom has also been discussed within models of the Nambu-Jona Lasinio class with the 
aid of vector interactions and Polyakov loops in Ref. || . 

The occurrence of inhomogeneous phases has also been studied in a number of different 
contexts including superconductivity (see for example Refs. [10, 11, 12, 13, 14, [H| |il| ), 



lower dimensional field theories and, in particular, the Gross-Neveu model with a large 
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number of fermions (see [17, 18, 19, 20, [2^, |2^]) that allows for exact integrability. 
Earlier studies regarding inhomogeneous phases were also carried out in the context of the 



Skyrme model [24, 25, 26, 27, [2^] and in four-fermion models p9|| . 

Strongly coupled theories are physically relevant also in applications of astrophysical 
and cosmological nature (neutron stars, black holes, cosmological phase transitions). In 
such cases it is natural to ask whether gravity may have some influence in the dynamics 
of phase transitions and, for this reason, the problem of generalizing strongly interacting 
fermion effective field theories to curved spacetimes triggered a large body of work (See 
Refs. |(| 0, H, fH H, |H H, gjj, H. Il> & M> HI HI B H for a Partial list of 



references). 

In curved space, mostly for reasons of technical nature, attention was limited to con- 
sider homogeneous spacetimes and condensates, i.e. spacetime curvature and condensate 
were assumed not to vary in space. Even in this relatively simple situation, the interesting 
conclusion was that gravity may, in fact, affect the critical temperature T c of the theory. 
Once the theory is immersed in a spacetime of positive constant curvature, in the mean field 
and large- N approximations, it was shown that, if temperature is fixed to some value T< 
smaller than the critical temperature T c , chiral symmetry is broken. Keeping the tempera- 
ture fixed to T< and increasing the value of the curvature, produces a phase transition into 
a chirally restored symmetry phase, indicating a clear analogy between thermodynamical 
and geometrical effects. 

Although this analogy is formally interesting, in concrete physical situations (for ex- 
ample in the case of a neutron star) the value of the curvature is small and geometrical 
effects only play a marginal role. In such a case, the system can be well described as if 
the spacetime were effectively flat, limiting the interest of the study of effective strongly 
interacting theories in curved space. However, this is not a universal situation. 

Clearly, geometrical effects become significant in the presence of strong gravitational 
sources. One natural example that comes to mind is that of black holes. This case, phys- 
ically very interesting, presents several complications absent in the cases of homogeneous 
backgrounds. The most important one is related to the fact that approximating the con- 
densate as a spatially constant function is simply not admissible. 

Targeting a solution to this problem, we have recently initiated the analysis of inho- 
mogeneous phases in curved spacetimes. In Ref. fltf , we developed a formalism that quite 
naturally allows to deal with this problem in regular ultrastatic spacetimes (the assump- 
tion of ultrastaticity was simply motivated by the fact that the special geometric structure 
allows a direct use of the imaginary time formalism necessary to include finite temperature 
effects). With some additional efforts, in Ref. [47], we have analysed the case of black holes 
exteriors and shown that chiral phase transitions occur outside the event horizon. More 
precisely, if the black hole is surrounded by a gas of strongly interacting particles, an inho- 
mogeneous condensate with a kink profile will form. Although the implementation was not 
straightforward, the idea is, in fact, simple. In a black hole geometry that is asymptotically 
flat, local quantities increase as the horizon is approached. If the asymptotic temperature, 
inversely proportional to the black hole mass, is lower than the critical temperature, then 
chiral symmetry is broken at infinity. According to the same logic, since the local tem- 
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perature diverges near the horizon, we may expect that, at some finite distance from the 
horizon, a transition to a chirally symmetric phase may occur. This suggests that for black 
holes of certain mass in equilibrium with a gas of strongly interacting particles, a bubble 
of high temperature restored chiral symmetry phase should surround the hole. For an 
evaporating black hole the same holds, since asymptotically the temperature vanishes as 
1/r 2 . This has been quantitatively demonstrated in [47] where self-consistent solutions for 
condensate have been constructed. 

In all cases analysed so far (for both homogeneous and inhomogeneous situations), the 
spacetime geometry was taken to be regular (in Ref. [46] we studied the case of regular 



ultr astatic manifolds, while in Ref. jlT]] we focused on the region outside the event horizon 
that is also singularity free). The goal of this paper is to analyse more deeply the connection 
between the appearance of inhomogeneous phases and curvature effects, in particular when 
the geometry presents singularities. Examples of this sort describe, for example, global 
monopoles spaces, geometries with wedges, string theory conifolds. Higher co-dimension 



brane models also present a similar structure [51]. 

The first step of the present work is to compute the effective action for a strongly 
coupled fermion effective field theory on a conical spacetime. Since the geometry presents 
singularities, a suitable regularization procedure has to be used to deal with this situation. 



One possibility is to use the approach developed by Cheeger in Ref. 4q| , where the spectral 
geometry of differential operators on singular spaces has been discussed. Although the 
problem can be discussed quite generally, here, for simplicity, let us specify the set-up and 
consider a quantum field \ propagating on a .D-dimensional spacetime ^# that presents a 
conical-type singularity. The effective action can be written in terms of the heat-kernel, J?T, 
of some differential operator defined on the spacetime manifold in question. Schematically, 
we may write 

jtr(t) ~ J2^/ 2ti/2 ~ D/2 + ^ lnt > ( L1 ) 

i 

where are the heat-kernel coefficients and the quantity SB is related to the value of the 
zeta function on the base manifold. In this case, the heat-kernel coefficients are understood 
as the principal part of integrals over the space manifold of geometric invariants and this 
is indicated by the hat. Using the above expansion, one may obtain the values C(0) and 
C'(0) related to the functional determinant of the operator in question. In general the zeta 
function may present poles and, as shown by Cheeger, by taking the principal part provides 
the appropriate treatment necessary in obtaining the asymptotic expansion of the heat- 
trace. In a sense, Cheeger's approach allows to deal directly with the singularity, however, 
in any practical implementation, truncating in the above expansion becomes necessary, and 
this affects the validity of the truncated expansion near the singularity. 

Another possibility, easier to implement and that we will use here, is to resolve the 
singularity by excising a small region around it and substituting the excised region with a 
regular cap (See Fig. ||) 1 . Adopting this way of regularization, allows to have a standard 



1 A similar type of regularization has also been used in the context of higher co-dimension brane models 
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expansion for J^(t), i.e. the coefficient 38 vanishes and the heat-kernel coefficients are the 
ordinary ones. 

Adopting this regularization sche- 
me, we will present the computation of 
the effective action for a strongly inter- 
acting fermion effective field theory on 
manifolds with conifold singularities in 
the next section. The method is anal- 
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Figure 1: Regularization by excision of the singular- 
ity. 



ogous to that presented in Ref. J46]. 
Once the explicit form of the effective 
action for the fermion condensate is ob- 
tained, we will proceed, in section || 
with the numerical implementation and 
construct self-consistent and spatially non-trivial solutions for the condensate. In partic- 
ular, we will see that even when the thermodynamical temperature is kept constant and 
below the critical value for the effective theory in flat space, chiral symmetry is always re- 
stored near the singularity. The condensate assumes a kink-like profile that vanishes near 
the singularity, thus separating the region of restored symmetry around the singularity 
from a region of broken symmetry. Our conclusions will be presented in the last section. 

2. Condensate Effective Action in Conifold-type Geometries 

In this section, we will consider the simplest tractable type of geometric singularity: that 
of a metric cone. Given a (d — l)-dimensional Riemannian manifold , a (metric) cone, 
< ^ / ^, is defined as the space M + x JV , In the following we will assume that the base of 
the cone, JV, is compact and smooth. The formal part of the analysis will deal with this 
general case, while in the subsequent numerical implementation we will specify the form of 
the base. In local hyperspherical coordinates the line element is 

ds 2 = g ij dx i dx j = dr 2 + r 2 d^ 2 , (2.1) 

with r being a radial coordinate, r £ R + . This is the geometry of a infinitely long cone 
and it is singular at the apex, r = 0, where a curvature singularity occurs. It is possible 
to cut the cone at some distance L, in which case the radial coordinate varies within the 
interval I = [0, L]. 

Finite temperature effects may be introduced using the standard imaginary time for- 
malism, thus considering the product of S 1 x ^€ jf where S 1 is a circle of radius 2ttT, and 
T = 1/(3 being the temperature in natural units. Anti-periodic boundary conditions will 
be imposed on the fermion fields on S . Generalization to scalars are almost trivial. 

We will model the strongly interacting fermion effective field theory with an action of 
the form 



S 



d d+1 x^\^V^ + ^(^) 2 



(2.2) 



fill |52| where the singularity was resolved by placing a brane at small distance from the singularity and by 
filling the interior by a regular cap. 
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where the spinor tjj has (4 x Nf x N c ) components, Nf and N c are the number of flavors and 
colors respectively (N = Nf xN c ), and A is the coupling constant. The dots represent terms 
with higher mass dimension. The matrices 7^ are the gamma matrices in curved space, 
and g = |Det<?y| with g^ given by ( |2.1| ). The theory is invariant under discrete chiral 
transformations, and possible phase transitions in relation to the breaking of the chiral 
symmetry can be discussed in terms of the appearance of a non vanishing condensate 
a = — ^(tptp): if the chiral symmetry is broken dynamically, a acquires a non-zero vacuum 
expectation value and a fermion mass term appears. The basics of the computation are 



the same as in our previous work Ref. [46], so, in recapitulating them, we will be brief. 

After bosonization, the effective action at finite temperature (per fermion degree of 
freedom), T, can be expressed in the large- ./V approximation as 

r = - J d d xj-g (^) + Tr In - a) , (2.3) 

where the determinant acts on spinor and coordinate space. Squaring the Dirac operator 
gives 



with 

6T = ^E Trln 



□ + | + <7 2 + e|<9<7| 



(2.5) 



where we have assumed the condensate to be symmetric around the axis of the cone, i.e. 
a = a(r). Zeta function regularization allows us to express 5T as 

<5r = 1 J d d x^g-(((0) ln£ 2 + C'(0)) , (2.6) 

where £ is a renormalization scale. The quantities £(0) and C'(0) are t ne analytically 
continued values to s = of the following function 

C ( s) = J_ y y f dtt s-i e -t^ Tie -K-A + ^ + e\ d a\) (2J) 

^ n e J 

In the above expression we have used the imaginary time formalism to introduce finite 
temperature effects (u n = 2ir(n + 1/2)//? are the Matsubara frequencies for the fermion 
fields), and we expressed the zeta function in terms of the Mellin transform of the heat- 
trace. 

In the following we will regularize the effective action by excising the singularity as 
this seems the most convenient way to deal with the subsequent numerical analysis. In 
this case the heat-kernel expansion will take the standard form and the radial coordinate 
is assumed to start from r = e, where boundary conditions will be imposed. Some algebra 
allows us to recast the above expression as 

C(a) = C+00 + C-(«) > (2-8) 
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where 



with 



C±(«) 



g(t) 



4 



dt 



t*- 1 e - tC7 ±je; e g(t) 

T(s) 



2Jitt 



git) 




(2.9) 
(2.10) 



Notice that in the regular part of the heat-trace we have factorized out the exponential 
according to the method described in Refs. fH^ , [50| ]. In the above expression s is regulating 
parameter assumed to lie in a region of the complex plane where all the above integrals are 
well defined. The standard expression for J^r eg (t) is given by 



( 4vr )rf/2 



■d/2 



(2.11) 



k=0 



where the coefficients are the heat-kernel coefficients associated with the differential 

operator in (2.3) and are integrals of polynomials in the curvatures on the regularized 
spacetime manifold and therefore are all regular. Coefficient with integer index refer to 
volume contributions while coefficients with semi- integer indices to boundary contributions. 



The first coefficient is 



(0) 



± 



1, while ^ 



0. Higher order coefficients will be given later 



as needed. 

Using the above expression for (±(s), we can express the effective action as 



r 



d A x 



V9[ 



( G 



^ l + r++r - 



(2.12) 



where T± are the contributions to the effective action coming from £±. For completeness, 
at the end of this section, we will illustrate how the computation can be performed also by 
using the full heat kernel expansion including the logarithmically singular term in (1.1). 

We will start assuming s to be in the region where the integrals are well defined, then 
analytically continue to s = and take the limit d = 3. The zeta function can be ex- 
pressed in terms of the following integral and its derivative with respect to the regularizing 
parameter s, 



J^(a, b, c) 



(47T 



dt 



j.s— a ^—ct—b/t 



(2.13) 



These can be computed explicitly, using the following expressions (see Ref. |55|1): 



J^(a, b, c) 

dJ r (a, b, c) 
ds 



(4vr)^r(s) 



r ( s ) (4vr 
dv 



d+l 

2 



(-l+o-s)/2 



K a . 



2Vbc 



(-l+a-a)/2 



(2Vte) (ln(0 +2^ (0) (s) 



(2y/bc\ 






v=a—s—l- 
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The above expression are valid under the following conditions: 



m > , 
Kc> , 

8fc(a-a)>l. (2.14) 

The first two requirements are trivially satisfied in our case. The third is instead dealt with 
in the usual way of zeta function regularization, i.e. by assuming that s lies in a region 
of the complex plane where the third inequality is satisfied and then proceed by analytical 
continuation. One may easily re-write the zeta function above as 

2 

e=± i=l 

where we have defined for notational convenience 

^\s) = fX /2 V(^±^,o,4) , (2-16) 

k=0 

= 2 ££(-i)^V(*±^,^ 1 ^ : ) , (2-17) 



Dividing now T± as the sum of the contribution coming from ((0), I± , and that 
coming from C'(O), I± , allows us to write 

r± = r2 ) + r£ ) . (2.18) 

A direct computation shows that only 3f±\s) gives a non vanishing contribution to C(0). 
We find 



r ± ) = 32^2 / *r 2 7det7 



ln£ 2 , (2.19) 



where det.js is the determinant of the metric on the base manifold JV . The above result is 
exact (no truncation in the heat-kernel expansion has been done and boundary contribu- 
tions also vanish). Also, notice that the ansatz used in the regular part of the heat-trace 

gives = (See Refs. jl^, |5Tj, ^6|). Therefore the second term in the above expression 

(2) 

vanishes. The coefficient ^ is given by 

tf[ 2) = —R. vpu m ivpcJ - —R UU R» V - —AR + -( ( a ± - ) <x" + a' 2 + -aa' 
x 18Q fiupcr lg0 p, v 12Q r ) r 

where the quantities R, R^ v and R^xp are the Ricci scalar and Ricci and Riemann tensors, 
respectively, for the geometry ( |2.1| ). All terms proportional to curvature invariants and not 
to the condensate will not, in fact, contribute to effective action for the condensate, in the 
sense that they will disappear in the equation of motion for a. The integrand in ( 2.19j ), 
i.e. the contribution to the Lagrangian density for the condensate, is regular in the limit 
r ->■ 0. 
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The term that requires more work is the contribution to the effective action that comes 
from the derivative of the zeta function, . Proceeding as described above, we find for 
the regularized effective action the following expression 



JL jdr r 2 Tdet^ j^i - Qc£ + *f ) In o\ 



+4^(-ir 



n=l 



4a 2 



^ 2 (n/3a ± )+^ 2) ^ (n/3cr±) 



+ 



(2.20) 



The dots represent higher order terms in the heat kernel expansions. They can be computed 
explicitly, but we won't report their form here. All higher order terms give contributions in 
the derivatives of the condensate of order equal or greater than four, so the expression above 
is sufficient when keeping the analysis to second order in the derivative of the condensate. 
Boundary terms can also be calculated easily and the first few terms are 



boundary 



00 

+*£(-!) 



(3/2) 
± 



n=l 



■ 2 V^^(l/2) (~ 1 



e -nl3a ± + 



ra/3 



(3/2) - n /3a ± 



+ ••■ 



where the dots represents, as before, higher order terms in the heat-kernel expansion. As 
before, one can easily verify the regularity of the expressions above. 

To conclude this section, we will show how the full heat-kernel expansion ([□]), includ- 
ing the logarithmic terms, can be used instead of the scheme adopted above. In this case, 
according to (1.1) the heat-kernel will consist of two terms. The polynomial part of the 
heat-kernel can be treated as above with the only difference being that the heat-kernel coef- 
ficients are understood as the principal part of integrals over the space manifold. Since the 
manifold has singularities, the integration in the heat-kernel coefficients, naturally, leads to 



divergences. As shown in Ref. [48], the relevant quantity for the heat-kernel asymptotics is 
the principal part of these diverging integrals that can be extracted in several ways. One 
possibility is to cut the range of integration at finite distance from the singularity and this is 
equivalent to the regularization by excision used above. The additional (logarithmic) term 



in (1.1) can also be included in the analysis, but some attention is necessary. The quantity 
is proportional to minus the residue of the zeta function at s = 0, — ResC^^O), that 
can, in turn, be related to the value of the zeta function evaluated on the base manifold 
analytically continued to s = —1/2. The logarithmic part of the heat-kernel expansion will 
lead to an additional term in zeta function that can be written as 



e=± i=3 



(2.21) 



- 8 - 



with 



iTi 3) (s) = -^^±(3/2,0,0) , (2.22) 

V47T 



^\s) = 2^ ^)(-l)»Jf(3/2, ^-,0) , (2.23) 

* n=l 

where 

JT (o, 6, c) = — ^ / *— In f . (2.24) 

(47r)~ J 1 l S J 

The analytical continuation of the logarithmic term can be done straightforwardly by using 
the following expressions 

^»= ^ r( - i +r ,) ^(- i+ — . 

(47r)~ 1 '""< 



2^2 



2 



d s (47r) ^ r( s ) 



(V 0) (-1 + a-s) - In ft) x 

(V ) (-1 + a - s) + V (0) (s) - In ft) + (-1 + a- s) 

whose analytical continuation can be performed as before. A straightforward computation 
shows that the contribution of the logarithmic term to £(0) vanishes, while the contribution 
to the effective action from C'(O) gives 



Jog _ _ ln2 



„ v7 e - In (2/3 2 )) J d d x^m ± . (2.25) 



(4tt)2/3 

3. Inhomogeneous Condensates in Monopole Geometries: Numerical Con- 
struction 

This section will be devoted to implement numerically the above formalism and to construct 
explicit solutions for the condensate. To have a concrete working model, we need to specify 
the geometry of the base. In order to make contact with a case of physical interest, we 
will choose the base manifold to be a sphere of non-unitary radius p. The curvature is 
R oc (1 - p 2 )p- 2 r- 2 , the space is non flat when p 7^ 1. In the following we will consider 



p < 1, therefore R > 0. In this case, the metric (2.1) describes the spatial section of global 



monopole [d3|. For this geometry, the heat-kernel analysis of the conformal Laplacian has 



been done, for instance, in Ref. [54]. 



The equation of motion for the condensate can be found by minimizing the effective 
action, that is by varying T with respect to the condensate a. In the previous section we 
have obtained the effective action to fourth order in the heat-kernel expansion {i.e. k = 4 
in ( [2.11D ), leading to a second order non- linear equation of motion for the condensate. This 
can be written, with some work, in the form of non-linear Schrodinger-like equation, whose 
explicit form is very lengthly and will be omitted here. 
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The structure of the solution may be anticipated by looking at the form of the ther- 
modynamic potential. At large distance from the singularity, where curvature effects are 
negligible, the thermodynamic potential U as (a) is expected to have a behaviour similar to 
flat space: below (above) some critical value of the temperature, chiral symmetry is broken 
(restored). This can be verified by letting a' = 0, taking the limit r — > oo of the effective 
potential, and by computing the thermodynamic potential by numerical integration of d a U 
with respect to a. The profile of the asymptotic potential is shown in Fig. |2| (left panel). 
The thermodynamic potential can also be computed locally and the result is illustrated in 
Fig. ^ (right panel) . The figure shows that even if chiral symmetry is broken asymptotically 
(right-most (purple) curve), chiral symmetry will gradually be restored as we move towards 
the singularity, i.e. the minima of the potential will gradually shift towards a configuration 
with vanishing a. The solution for the condensate can be found using standard numerical 
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Figure 2: The figure on the left illustrates how the local thermodynamic potential changes as the 
temperature increases. The right-most (left-most) curve shows the potential for temperature smaller 
(higher) than the critical one. The figure on the right illustrates how the local thermodynamic 
potential changes as the singularity is approached. The right-most curve shows the potential at 
large distance from the singularity for a set of parameters that correspond to a chirally broken 
symmetry phase (The minima of the potential is non vanishing). The left- most curve shows the 
potential close to the singularity (r = 0.1) for which the potential has a vanishing minima at a = 0, 
leading to a chirally restored symmetry phase. For both figures we set I — 10 6 and A = 10 -2 . 



methods, and here we have used a fourth order Runge-Kutta one. In solving the equation 
for the condensate numerically, there are several things that require some care. One has 
to deal with the infinite summations over Bessel functions with argument proportional to 
the condensate. The present situation differs from the black hole case studied earlier [47] 
as described in the following. At large distance curvature effects are negligble and the 
argument of the Bessel function is not small. In this case, one may conveniently truncate 
the summation due to the fact that the Bessel function decay exponentially for large ar- 
guments. Near the singularity the contribution from the curvature is, instead, large, and 
the argument of the Bessel functions is again not small. In the intermediate region for 
r sufficiently large, we require that, when the argument of the Bessel functions becomes 
small, the summands are expanded for small arguments and full resummation over n is 
performed. We then proceed by matching the solutions in the truncated and resummed 
domains. With the equations of motion for the condensate explicitly written, we solve 
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Condensate profile found 
minimizing the effective 



Figure 3 

numerically 

tion (2.3), for four indicative 
base manifold radius (Left to 



of 



by 

ac- 
the 



them by requiring regularity for the solution. The boundary conditions on the solution are 
set numerically by requiring, near the singularity, that the condensate is at a minima of 
the potential and by fine tuning the value of the derivative to match the asymptotic value 
of the solution with the minima of the potential at infinity. 

The results for the condensate pro- 
file are illustrated in Fig. || for sample 
values of the parameters showing that 
near the singularity the condensate van- 
ishes and chiral symmetry is restored, 
while, as we gradually move away from 
the singularity, the condensate assumes 
a non vanishing expectation value and 
chiral symmetry breaks. Notice that con- 
trary to the black hole case we have 
studied earlier p7| , in the present case, 
the thermodynamical temperature is con- 
stant uniformly along the spatial sec- 
tion of the space and the effect of break- 
ing/restoration of the chiral symmetry 
is solely due to genuine curvature ef- 
fects. As in the black hole case fl47|| , 
this indicates that a bubble of chirally 
restored symetry phase surrounds the 
singularity. 

The analysis can be extended to higher 
orders by including higher order terms 
in the heat-kernel expansion. The algebra becomes rather cumbersome, but it is possible 
to handle, with some work, all the computations automatically by using any computer 
algebra program. After including terms up to fourth order derivative of the condensate, 
we implement our analysis by perturbing the background solution, a = a^g + 5a, and 
expanding the higher order equation for a keeping terms up to second derivative of the 
perturbation. Proceeding in this way, we are assuming that the corrections due to higher 
order terms only produce small changes in the background solution. This assumption will 
be verified a posteriori after the solution is obtained. Once higher order derivatives of the 
perturbation are dropped out, the equation for the 5a becomes again a second order one 
with source term and can be solved by standard methods. This approach is similar to the 



values 

right: p = 
0.65 (Green), 0.57 (Blue), 0.53 (Black), 0.48 (Red)). 
We set to £ = 10 6 , A = 10~ 2 . The superposed panel 
shows the solution for the condensate when higher 
order terms included superposed over the background 
solution. The continuous line refers to the back- 
ground solution corresponding to p — 0.65 (Green) 
and the dotted (red) solution corresponds to the 
higher order one. 



case of black holes we have addressed earlier [47] and, as in that case, the corrrection to 
the solution only produces very small distortions that decrease with the thickness of the 
kink. A sample plot is shown in the small panel in Fig. H. 



4. Conclusions 



In this paper we have continued the analysis of the chiral symmetry breaking of strongly 
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coupled fermion effective field theory in curved space, extending our previous work to 
include conifold-type geometries. 

The results, in agreement with our intuition, show that, near the singularity, the 
condensate vanishes and chiral symmetry is restored. As we gradually move away from the 
singularity, the condensate assumes a non vanishing expectation value and chiral symmetry 
breaks. This indicates that the singularity will be surrounded by a bubble of chirally 
restored symmetry phase, in much the same way as black holes (See Ref. ^7|). We have 
illustrated this explicitly in the tractable example of a conifold geometry with generic base 
manifold (the generalized metric cone). More precisely, we have considered a strongly 
interacting fermion effective field theory propagating on the conifold and constructed an 
explicit solution for the condensate. The formal analysis has been carried out in general, 
and we have evaluated the effective action for the condensate introducing finite temperature 
effects by means of the standard Matsubara formalism. This set-up allowed us to keep the 
genuine thermodynamical temperature constant and to clearly single out any non-trivial 
effect of the curvature. The evaluation of the effective action was carried out using a 
sophistication of the method we have described in Ref. [46] and that makes use of zeta 
function regularization. Due to the presence of the singularity the regularization requires 
some care, however zeta function regularization proves to be very adequate in this case. 

Several generalizations may be considered. First of all, one expects that the present 
situation is more general. That is in the vicinity of any spacetime singularity the same 
occurs. This can in principle be treated by studying the heat-kernel asymptotics in space- 
times with different singular structure. Although computational complications may arise, 
the procedure presented in this paper that regularizes the geometry by excising a small 
region around the singularity should work with little modifications. Generalizations to 
more sophisticated models may also be possible with some effort, the most interesting ones 
being the inclusion of gauge degrees of freedom and of a chemical potential. 

It was our goal to understand how condensates and chiral symmetry may be affected by 
curvature singularities and for this reason we considered a purely geometrical background 
(2.1). However, the results may be easily extended to cases where a non-trivial radial de- 
pendence of the lapse function occurs, and this can be handled by using the same conformal 
techniques that we have adopted to study the case of black holes. In this case, if the local 
temperature increases as the singularity is approached we expect the same phenomena of 
chiral restoration described here to occur. 
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